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ABSTRACT 

An elementary derivation is given for the "Peierles substitution" used in projecting dy- 
namics in a strong magnetic field onto the lowest Landau level. The projection of wavefunc- 
tions and the ordering prescription for the projected Hamiltonian is explained. 
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The "Peierles substitution" 1 advances calculations for the following problem. Consider a 
charged particle with mass m moving on the plane in a constant magnetic field B perpendic- 
ular to the plane, and also subject to planar forces arising from the potential V(x, y). Motion 
along the z-axis (along B) is ignored and it is assumed that B is sufficiently strong and V 
sufficiently weak so that the lowest Landau level retains its identity in the presence of V. The 
question one wants to answer is how the energy-degenerate states in the lowest Landau level 
are modified by the interaction with V. The result in the approximation of strong B and 
weak V is that the energy eigenvalues become 

where B/2m is the lowest Landau level energy [we set h, c and particle charge to unity] while 
e n are eigenvalues of the operator obtained from V(x,y) by the Peierles substitution, 

V(p,q)\n) = e n \n) (2) 

where p and q are canonically conjugate (with 1/B playing the role of Planck's constant). 

M = ^ (3) 

The Peierles substitution has been used for over half a century; a recent application is 
to the Azbel ,2 -Hofstadter 3 problem where V(x, y) = V cosx + V cosy. But its justification 4 
remains a "most difficult step." 3 

In this Letter we draw on our experience with Chern-Simons quantum mechanics 5 to give 
an immediate derivation of the Peierles substitution, and to answer two immediate questions: 
first, how do the wavefunctions ip(x,y) which depend on two commuting coordinates x and 
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y reduce to wavefunctions depending only on the single coordiante ql Second, how does one 
resolve the operator ordering ambiguity inherent in (2) when the potential V(x,y) becomes a 
function of the non-commuting operators p and ql We answer these questions and explicitly 
verify our prescription when V is rotationally symmetric: V = V(x 2 + y 2 ). 

Our derivation relies on the following fact about quantizing a system governed by a 
Lagrangian that is linear in time derivatives. Let £ l (i = 1, 2) be a two-component quantity 
with Lagrangian 

/- §£eV ViO (4) 

ij 

Here g is a constant, the anti-symmetric tensor = ( ^ J J , and V a function 

V / ij 

depending on but not on its time-derivative £\ Although lacking the usual kinetic term, 
quadratic in time derivatives, (4) gives rise to non-singular equations of motion, which may 
be quantized in a Hamiltonian framework, provided the fundamental quantum commutator 
is taken as 

[e,£ J ] =~e ij • (5) 
9 

While (5) is a standard result 6 about systems like (4), it is not well-known, hence in the 
Appendix we sketch a derivation. 

Consider now the magnetic problem described by the Lagrangian 



L m = ^mr 2 + | r x r - V(r) (6) 



Dynamics is confined to the plane, r = (x, y), the ^-coordinate is suppressed and the vector 
potential A giving rise to the constant magnetic field B = V x A is taken in the rotationally 



symmetric gauge A l (x) = —e^r^. 



The Hamiltonian corresponding to (6) is the familiar 



expression 



#m(p,r) 



1 



(p-A(r)) 2 + y(r) 



(7) 
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In the absence of V, the spectrum is well-known: the infinitely degenerate lowest Landau 
level has energy B /2m and the higher levels are separated by gaps of magnitude proportional 
to B/m. When V enters weakly and B is strong, the degeneracy is lifted, but the pattern 
remains unchanged, and the lowest level may be isolated by setting m to zero, provided the 
energy eigenvalue B /2m is subtracted ( "renormalized" ) . Thus the structure of the lowest level 
is governed by the m = limit. 

[The analogy to gauge theories with Chern-Simons terms is now evident: the kinetic and 
potential terms in (6) are point-particle analogues of the conventional Maxwell/ Yang-Mills 
Lagrangian; the second, magnetic term models the field theoretic Chern-Simons interaction. 
Setting m to zero in the quantum mechanics problem parallels the passage from a topologically 
massive gauge field theory (with the Maxwell/ Yang-Mills and Chern-Simons terms) to a pure 
Chern-Simons field theory. Indeed Ref. [5] is devoted to an illustration of these field-theoretic 
issues in a quantum mechanical context. But here now we pursue the actual physical relevance 
of the quantum mechanical model.] 

The m — > limit may be discussed in the Lagrangian (6) or Hamiltonian (7) framework. 
In view of the previously stated result (4) and (5), the Lagrangian discussion directly leads 
to the conclusion that the symplectic structure of the reduced Lagrangian 



— r x r — V(r) 



(8) 



enforces the commutation relation 





(9) 
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and the reduced Hamiltonian H , which is obtained from L by the usual Legendre transform 

H = ^f-r-L = V(r) , (10) 

becomes evaluated with non-commuting arguments in view of (9). An alternative, more 
complicated, derivation starts with H m of (7): the limit m = is achieved by imposing 
quantum mechanically the constraint 

= p - A = mr (11) 

H reduces to H and (9) emerges as a consequence of constrained quantization. 5 ' 6 

We have thus arrived at a derivation of the Peierles substitution: the perturbed structure 
of the lowest Landau level is governed by the effective Hamiltonian 

H = V(p,q) (12) 

where the non-commutativity of the arguments follows from (9) (expressed as in (3)). 
When V = V(r 2 ) = V(p 2 + q 2 ) it is temptng to write the eigenvalues of H as 

en~v(l(n+]\) (13) 



where (n + ^) are the harmonic oscillator eigenvalues. However, this naive substitution 
ignores operator ordering issues, which produce "corrections" to (13). 

To study the operator ordering we first describe the behavior of wavefunctions when the 
limit m — > is taken. Note that in the limit, phase space is reduced from four dimensions 
{Px,P y ,x,y) to two (p, q); this is seen also as a consequence of the constraint (3). Wave 
functions depend on half the phase-space variables: two before the reduction and one after, 



but if normalization is maintained, one argument of ip m (x, y) cannot simply disappear when 
m vanishes. Detailed analysis 5 gives the following story. In quantizing the reduced theory, 
the holomorphic polarization is chosen for wavefunctions: 7 the single combination of the two 
phase space variables, on which wavefunctions are taken to depend, is the non-Hermitian 



combination (x + iy). Thus we form the operators 

a = \j^{x-iy) , a! = yj^{x + iy) , [a,a f ] = l (14) 
Coherent states provide a basis 

(a\af = (a\a (15) 
and states are described by wavefunctions that depend on a. 

HV> = V(«) (16) 

The operator a) acts on these functions by multiplication, a by differentiation; the adjoint 

relationship between the two is maintained by virtue of a non-trivial measure. 

{a\a)\^ = aip(a) 



d 

(a\a\ip) = ^V>(a) 

— / da*dae~^ \a)(a\ = I 
2m J 1 A 1 

da* da = — dxdy . 

2ni 2tx y 



(17) 



(18) 



One is also interested in number states 

a^a\n) = n\n) (19a) 
which are described within the holomorphic representation by 



(a\n) = -= . (19b) 
ynl 
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The phase-space reductive m — > limit onto the lowest Landau level, when applied 
to wavefunctions, works in the following manner. The (x, y)-dependence in the fmite-m 
eigenf unctions of H m for the lowest Landau level, are presented in terms of complex variables 
z = ^/Bj2(x + iy), z* = ^/B]2{x - ty) 

tp m = ^m(z,z*) (20a) 

When m is set to zero, the limiting form becomes 

/ d \ 1/2 

limj m (z,z*)=( — ) e-il 2 ! 2 ^) (20b) 

with ip being a holomorphic eigenfunction of H . Note that the probability densities transform 
properly so that the norm is maintained. 

d 2 r |^ m (r)| 2 — f dxdye-^mz)\ 2 = ^! e -M 2 \^ a) f (21) 

m— *-0 /7T Z7TZ 

(All this may be easily and explicitly checked for the solvable case when V is a harmonic 
oscillator potential. 5 ) 

The ordering prescription for Hq = V must be such that the matrix elements e = (n\V \n) 
of V, when computed before the phase-space reduction as integrals over the function V(x, y), 
coincide with the evaluation after reduction, when V becomes an operator between the number 
states | to) in (19a). This prescribes that the operator H be expressed as the anti-normal or- 
dered (i.e. all a's to the left, a^'s to the right) form of V when V(x, y) is written as V(z, z*) and 
z* is replaced by a and z by aL (To see the need for anti- normal ordering, note that the coher- 
ent state resolution of the identity in (18) implies that a k (a^Y = J e~^ z ^\z) (z*) k z e (z\.) 

7 



To verify this explicitly, we consider a potential V = V(r 2 ) as a function of r = p +q ~ 
-^a^a, a potential where the tilde indicates that the ordering is yet to be performed. Then 
mH m may be written as 

mH m = i(p- A) 2 + mV (22) 

and the energies e n of (1) can be evaluated by lowest order perturbation theory in m. The 
ground state of the unperturbed Hamiltonian (p — A) 2 /2 carries energy B/2; it is infinitely 
degenerate with wavefunctions 



"" (r) = v?UJ • n = ' 1 "" (23) 

Degenerate perturbation theory must be employed, but fortunately a rotationally symmetric 
V is already diagonal in the representation (23). Thus we find 

(jj \ l+n poo 
| J J dr 2 r 2n V(r 2 )e-^ 2 (24) 



If V is represented as 



it follows that 



V{r 2 ) = / dae~ ar v(a) (25) 



2 x " 1_n 



e n = j da{l + -a v(a) (26) 

On the other hand, in the reduced problem we have the operator Hamiltonian 

V ~ y rfae-^^^fa) (28a) 

To anti-normal order, we expand the exponential 

"=/*»wEf(-i)"(W 

= Jdav(a) (l+ jr ("^-) ( atfl + iV ) ( ata + ^ " • • • ( ata + 2 ) ( ata + 

(28b) 



The number states \n) are eigenstates with eigenvalue 



/ 



N=0 v 7 



N 



(N + n)\ 



(29) 



Summing the series reproduces (26). These expressions, (26) and (29), encapsulate all the 
"corrections" to the naive replacement in (13). 



Dimensional reduction of phase space in the presence of strong magnetic fields has re- 
cently also been studied by Iso, Karabali, Sakita, Sheng and Su, 8 as well as by Levit and 
Siran. 9 We thank some of these authors for discussion. Levit and Sivan give a semi-classical 
(WKB) analysis of ordering corrections to Hq; their prescription in general disagrees with 
ours, agreement exists only for the first two terms of a large-n expansion. We are grateful 
to C. Callan and D. Freed for acquainting us with the Peierles substitution and M. Berry for 
suggesting the perturbative calculation. 



After completing this work, we searched the literature and discovered that most of our 
points were made almost a decade ago by S. Girvin and T. Jach, Phys. Rev. D 29, 5617 
(1984); see also S. Kivelson, C. Kallin, D. Arovas and J. Schrieffer, Phys. Rev. B 36, 1620 
(1987). 
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APPENDIX 

We consider the following Lagrangian that governs dynamics for a multiplet of variables 
£ l and is linear in the velocities £*. 

L = a i (Z)?-V(Z) (A.l) 

Here ai and V are functions of but not of £ l , and the summation convention is used. L of 
(A.l) is a generalization of (4), where aj(£) = f e y'C" ? - Although (A.l) lacks the usual kinetic 
term, quadratic in velocities, it is in fact the generic form for a Lagrangian. The point is that 
familiar Lagrangians like 

L = \mq 2 - V(q) (A.2a) 

may be expressed as 

2 

L(p,q)=pq-H(p,q) , H(p, q) = |- + V (q) (A.2b) 

and the Euler-Lagrange equations for L(p,q), with p and g taken as Lagrangian variables, 
reproduce the equations of motion of (A. 2a). But (A. 2b), is of the form (A.l) with £ l = (^). 
The equations of motion that follow from (A.l) read 

f'Me = %^- (a.3) 

^(0 = 1^(0- (A.4) 

Next, we assume that the matrix possesses an inverse f lj — if it does not, the subsequent 
development is more complicated, 6 but unnecessary in our case (4). Then (A.3) implies 
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The Hamiltonian for (A.l) is constructed by the usual Legendre transform 

r)T 

H = ^C-L = V(0 (A.6) 

Hence if the aim is to reproduce (A. 5) by bracketing with the Hamiltonian, it must be that 
the bracket isnon-vanishing since H = V(£) depends only on £*. Evidently, 

i[y(o,c] = ^^-i[e,e} = e (a.?) 

and (A. 5) is regained provided 

e-e\=if ij (0 (A.8) 

This gives the derivation of (5). 
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